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Abstract 

In multiply-charged ion plasmas, a significant number of electrons may occupy high-energy orbitals. 
These “Rydberg” electrons, when they act as spectators, are responsible for a number of satellites of 
X-ray absorption or emission lines, yielding a broadening of the red wing of the resonance lines. The 
contribution of such satellite lines may be important, because of the high degeneracy of the relevant 
excited configurations which give these large Boltzmann weights. However, it is in general difficult 
to take these configurations into account since they are likely to give rise to a large number of lines. 
We propose to model the perturbation induced by the spectators in a way similar to the Partially- 
Resolved-Transition-Array approach recently published by C. Iglesias. It consists in a partial detailed- 
line-accounting calculation in which the effect of the Rydberg spectators is included through a shift and 
width, expressed in terms of the canonical partition functions, which are key-ingredients of the Super- 
Transition-Arrays model. The resulting method can a priori be used in any detailed-configuration/line- 
accounting opacity code. 


Keywords 

hot plasma - atomic physics - highly excited states - satellite lines 

PACS 

32.70.Cs, 32.70.-n, 32.80.Ee, 32.80.Zb 


1 Introduction 

Highly charged ions may have electrons in their high-lying loosely-bound orbitals. These “Rydberg” 
electrons, when they act as spectators, yield satellite lines that may not be clearly separated from the 
resonance line: they contribute as a broadening of its red wing. Some satellite lines are sufficiently 
strong that they can be used for diagnostics purposes. In ICF (Inertial Confinement Fusion) plasmas, 
electron temperature T and density Ae of the fuel are inferred from the X-ray emission of H-like and 
/ or He-like lines of impurity ions, e.g., Ar ions m), which are enhanced and broadened by unresolved 
satellite lines. Concerning the Hohlraum emission, satellites of 3d-4f transitions of Ni-like gold ions form 
a quasi-continuum in the 2-4 keV energy range [5]. In another context, a large number of satellite lines 
due to 4d-4f transitions contribute significantly to the EUV emission of Xe over the range [10-17] nm [5] 
(see Fig. [^, which is important for lithographic applications. 

The contribution of dielectronic satellite transition arrays, with spectator electrons in shell n=5, was 
studied by Bauche-Arnoult et al. for the case of a tantalum plasma |3]. The main mechanism for the 
formation of doubly-excited states is the dielectronic recombination. However, for high electron densities, 
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Figure 1: (Color online) Effect of passive subshell nd^, n=5 to 9, on transition array [Kr] 4d® —>■ 4d®4f^. 
The calculation was performed with Cowan’s code (RCN/RCN2/RCG). 


those excited states can also be populated by successive collisional and radiative excitations from the inner 
shells. In general, the study of the corresponding satellite lines requires a non-LTE (local thermodynamic 
equilibrium) modeling, in which all microscopic processes are accounted for in a master equation in order 
to determine the population of states. Most of the doubly-excited states involved lie above the ionization 
potential of the next ion, so that they may broaden some transition arrays by autoionization. Biedermann 
et al. measured the X-ray emission from high-n spectator electrons following radiative stabilization 
in dielectronically excited to Ba'^^^ ions produced by an electron beam ion trap. They observed 

an enhanced emission due to transitions from n >3 to n=3. 

In LTE, the contribution of those satellites may be important because of the large Boltzmann sta¬ 
tistical weight of the relevant excited configurations (due to high electron temperature and/or high 
degeneracy). In fact, even if their individual probabilities are rather small, a large number of excited 
states may fill gaps in the spectrum [^. In general, however, it is difficult to obtain a detailed description 
of those satellites because the number of electron configurations may become huge and they can give rise 
to a large number of lines (especially when open subshells with large angular momentum i are involved). 
In that case, detailed-line-accounting (DLA) calculations become intractable. 

A short review of semi-empirical and approximate models for treating high-n spectators is presented 
in section!^ These methods were proposed to address the aforementioned issue in a pragmatic manner. 
In section^ the detailed opacity code SCO-RCG mil] is briefly described. The statistical approach 
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that is currently used in the code to model those Rydberg spectators is explained. In section]^ we show 
the limit of this approach by studying the F-like 2p—)■ 4d transition of an iron plasma at r=182 eV and 
Nfi = 3.1 10^^ cm“^. These conditions correspond to a recent experiment performed by J. Bailey on the 
Sandia Z machine [3]. In section]^ we propose a new approach which allows one to improve the modeling 
of Rydberg spectators in the opacity code. 

2 Semi-empirical approaches and approximate models 

In Detailed-Configuration-Accounting codes, the Rydberg electrons can be only taken into account par¬ 
tially or approximately, because they are potentially responsible for a dramatic increase of the number of 
electron configurations. Thus, it is in general difficult to take those satellites into account, as they give 
rise to a large number of lines. In previous applications these contributions are either partially discarded, 
or described semi-empirically. Over the years, semi-empirical approaches and approximate models were 
proposed to solve that problem. In the middle of the eighties, Busquet et al. introduced a broad empirical 
red wing nniiii] in order to reproduce the M-shell gold spectra in the presence of myriad unresolvable 
satellite lines. At the same period, Goldberg and Rozsnyai m proposed to add a Gaussian width to the 
lines, based on statistical fluctuations of the populations of spectator electrons in high-lying states. Such 
an approach does not reproduce the red shift but accounts approximately for the way Rydberg spectators 
fill gaps between clusters of lines [13] . The correction to the variance reads (see Appendix): 

^dielectronic — ^ ^ ^ j ’ (^) 

where Vij denotes the mean electrostatic interaction energy between subshells i and j, Np the average 
electron population of subshell Upip and gp its degeneracy. Several years later, Honda et al. jTH [T^ ITS] 
performed intensive MCDF (Multi-Gonfiguration Dirac-Fock) calculations (using the GRASP2 code (TT] 
fTS] l and proposed to extrapolate the satellite emission spectra of gold ions from a single spectator to 
an arbitrary number of spectators. Such an approach is interesting, but inadequate when two or more 
spectator electrons are in the same shell. 

More recently, Iglesias et al. presented the PRTA (Partially Resolved Transition Array) model [T31 [301 
lU, which consists in representing a transition array by a small-scale detailed calculation that preserves 
the known properties of the transition array (energy and variance [33] [33] |3T]) and yields improved 
higher-order moments. In the PRTA approach, open subshells are split in two groups. The main group 
includes the active electrons and the electrons that are strongly coupled to them. The other subshells 
are relegated to the secondary group. A small-scale DLA calculation is performed for the main group 
(omitting therefore the subshells in the secondary group) and a statistical approach for the secondary 
group assigns the missing variance to the lines, which can be obtained analytically |23j . For instance, 
one may replace the transition array C C 

f C' = K2L®(3s)2(3p)2(3d)3(4s)i 
\ C" = K2L®(3s)2(3p)H3d)3(4s)i(5s)i 

which contains 102675 lines, by the “reduced” transition array C ^ C' 

f C = K2L®(3s)2(3p)2(3d)3 
I C" = K2L8(3s)2(3p)i(3d)3(5s)i 

which has only 26903 lines. The missing variance, i.e., the contribution of 4s^) is added as an additional 
broadening to the lines of C —> C'. The additivity of the contributions of active and passive subshells to 
the total variance is valid under the assumption that the radial integrals (Slater and spin-orbit) are the 
same for the initial C and final C" configurations. We can see for the example in Fig. [^for the transition 
3p3/2 5s of Hg, that the simplified calculation is very close to the exact one, although the transition 
array is very asymmetrical. 


3 



Figure 2: (Color online) Comparison between two SCO-RCG calculations relying respectively on DLA 
and PRTA treatments of lines for transition array 3p3/2 —>■ 5s in a Hg plasma at T=600 eV and /9=0.01 
g/cvo?. The DLA calculation contains 102 675 lines and the PRTA 26 903 lines. 
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3 The SCO-RCG code 


When lines coalesce into broad unresolved patterns due to physical broadening mechanisms (Stark effect, 
electron collisions, autoionization, etc.), they can be handled by global methods [25]. On the other hand, 
some transition arrays exhibit a small number of lines that must be taken into account individually. Those 
lines are important for plasma diagnostics, interpretation of spectroscopy experiments and for calculating 
the Rosseland mean, which is very sensitive to the hollows in the spectrum. For that reason, we have 
been developping over the past five years the “hybrid” opacity code SCO-RCG jam, which combines 
fine-structure calculations with statistical methods, assuming local thermodynamic equilibrium. 

3.1 Structure of the code 

Data required for the calculation of lines (Slater, spin-orbit and dipolar integrals) are provided by SCO 
(Super-configuration Code for Opacity) |26j . which takes into account plasma screening and density 
effects on the wave-functions. Then, level energies and lines are calculated by an adapted routine (RCG) of 
Gowan’s atomic structure code m performing the diagonalization of the Hamiltonian matrix. Transition 
arrays for which a detailed treatment is not required or impossible are described statistically, by UTA 
(Unresolved Transition Arrays) |22|, SOSA (Spin-Orbit Split Arrays) [24] or STA (Super Transition 
Arrays) [5S] formalisms implemented in SCO. 

Some other approaches were proposed to overcome the limits of statistical methods. For instance, in 
order to improve the accuracy of the STA approach, Hazak and Kurzweil |29j developed an interesting 
method, named CRSTA (Configurationally Resolved Super Transition Arrays); however, its degree of 
refinement is limited to the configurations, the terms structure being included statistically through UTA 
widths and shifts [3(1] . 

3.2 Probabilities of levels, configurations and super-configurations in the hy¬ 
brid formalism 

The total opacity is the sum of photo-ionization, inverse Bremsstrahlung and Thomson scattering spectra 
calculated by SCO code and a photo-excitation spectrum arising from contributions of SCO and Cowan’s 
codes in the form 


K {hv) 


1 M TTe^h 
47reo A me 


fx^x'Vx^x^x'ihv), 

X^X' 


( 2 ) 


where h is Planck’s constant, M the Avogadro number, eg the vacuum polarizability, m the electron 
mass, A the atomic number and c the speed of light. 7^ is a probability, / an oscillator strength, 4>(/ii/) 
a profile and the sum A —>■ A' runs over lines, UTAs, SOSAs or STAs of all ion charge states present 
in the plasma. Special care is taken to calculate the probability of A, which can be either a level aJ 
(a is a set of quantum numbers that uniquely specify the level and the coupling scheme used to obtain 
the total angular momentum J), a configuration C or a super-configuration S, because it can be the 
starting point for different transitions (DLA, UTA, SOSA or STA). In order to ensure the normalization 
of probabilities, we introduce three disjoint ensembles: V (detailed levels aJ), C (configurations C too 
complex to be detailed) and S (super-configurations S that do not reduce to ordinary configurations). 
The total partition function then reads 


Utot = U {V) + U (C) + U (5) with D n C n 5 = 0, 


( 3 ) 


where each term is a trace over quantum states of the form Tr 


„-/3(A-/.JV) 


H being the Hamiltonian, 


N the number operator, /x the chemical potential and /? = 1/ (ksT). The Hamiltonian takes into account 
the kinetic energy and the bound-bound, bound-free and free-free electron-electron interaction energies. 
The probabilities of the different species of the A-electron ion are 


5 





1 


(2J + 


( 4 ) 


for a level belonging to V, 


VaJ = 


Utot 


Vc 


1 


^ (2J+ 
ajec 


for a configuration that can be detailed, 


( 5 ) 


Vc = j^gc 

Otot 

for a configuration that can not be detailed {i.e. belonging to C) and 
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for a super-configuration. In the case where the transition C —?> C" is a UTA that can be replaced by a 
PRTA (see Fig. [^, its contribution to the opacity is modified according to 


fc^C Vc '^C^C'ihv) « ^ faJ^a'J' V^J aj^a'J' ( 8 ) 

where the sum runs over PRTA lines aJ —>■ a'J' between pseudo-levels of the reduced configurations, 
faJ^a'J' is the corresponding oscillator strength and j, is the line profile augmented with the 

statistical width due to the other (non included) spectator subshells. The probability of the pseudo-level 
aJ of configuration C reads 


with 


Vf^j — 


{2J + 1) e’ 


f' V“Q: j 




X V 


Y] VaJ = Vc, 

aJGC 

where Vc is the probability of the genuine configuration given in Eq. Q. 


(9) 


( 10 ) 


3.3 Interpretation of experiments 

The SCO-RCG code has been successfully compared to spectra measured in experiments performed on 
laser-produced plasmas, for example, (Fig. and a Z-pinch, for example, (Fig. [Mj facility. 

As mentioned in Sec. the quantity which is measured experimentally is the transmission, related to the 
opacity by Beer-Lambert-Bouguer’s law: T{hv) = ^-pLK{hu) ^ where L is the thickness of the sample. This 
relation between transmission and opacity is valid under the assumption that the material is optically 
thin and that re-absorption processes can be neglected. 

In the present version of the code, configuration mixing is limited to electrostatic interaction between 
relativistic sub-configurations {n£j orbitals) belonging to a non-relativistic configuration {ni orbitals), 
namely “relativistic configuration interaction”. That effect has a strong impact on the ratio of the two 
relativistic substructures of the 2p—>• 3d transition, visible around A=13 A in Fig. 
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Figure 3: (Color online) Interpretation with SCO-RCG code of the copper spectrum (2p—3d transitions) 
measured by Loisel et al. EIIIMIISS]. The temperature is T=16 eV, the density p=5 10 ^ g/cm^, the 
areal mass pL=AQ pg/cw? and the instrumental width 2.8 eV. The gray set of lines corresponds to the 
calculated spectrum before accounting for the intrumental broadening. 



Figure 4: (Color online) Interpretation with SCO-RCG code of the iron spectrum (2p—)■ 3d transitions) 
measured by Bailey et al. [34]. The temperature is r=150 eV, the density p=0.058 g/cm^, the areal 
mass pL=32 pg/cm^ and the instrumental width 1.8 eV. 
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Figure 5: (Color online) Effect of passive subshell nd^, n=4,5,6,7,8 and 15, on transition array [Ar] 3d® 
—>■ 3d®4f^. The calculation was performed with Cowan’s code (RCN/RCN2/RCG). 

3.4 Previous approach used in SCO-RCG 

Figures[5]and[^illustrate the effect of adding a passive subshell nd^ and nR respectively, with n=4,5,6,7,8 
and 15, to the transition array 3d® —)■ 3d®4R. The calculation was performed with Cowan’s suite of codes 
(RCN/RCN2/RCG). We can see that when the spectator electron is 4d, the spectrum of 3d®4d^—>■ 
3d®4d^4R is very different from the one of 3d®—3d®4R. This is due to the fact that the spectator 
electron 4d is strongly coupled to the active ones 3d and 4f. When the electron is pushed away (up to 
15d and 15f), the spectrum becomes closer and closer to the resonance lines. 

In SCO-RGG [7l[8j, orbitals are treated individually up to a certain limit beyond which all remaining 
orbitals (up to 15t) are gathered in a single super-shell a, called attributively “Rydberg super-shell” in 
the following. The grouped orbitals are chosen so that they weakly interact with inner orbitals. The last 
orbital of the super-shell corresponds to the one with the highest n and £ values for which a wavefunction 
can be calculated within the ion-sphere model. This corresponds to the pressure-ionization limit beyond 
which orbitals are assumed to be in the continuum. The grouped orbitals can also be defined by the user. 
In the latter case, it is useful to have an estimate of the last populated shell. For high values of n, the 
members of the series are found to be closer and closer together. When the finite width of the lines is also 
taken into account, then, for the highest levels, the spacing between two neighboring lines can be smaller 
than the line’s widths. Beyond this limit the lines cannot be resolved as separate lines any longer, but 
rather generate a smooth continuum. For the hydrogen-like Lyman series, the highest line that is still 
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Figure 6: (Color online) Effect of passive subshell nf^, n=4,5,6,7,8 and 15, on transition array [Ar] 3d® 
—> 3d®4f^. The calculation was performed with Cowan’s code (RCN/RCN2/RCG). 
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resolvable was calculated by Inglis and Teller [3^ assuming a Stark broadening. The expression of the 
maximal principal quantum number nmax, determined by Inglis and Teller from a simplified screening 
model in a very dilute hydrogen plasma, was extended to non-hydrogenic system for instance by Alaterre 

m- 


logio = 23.26 + 4 logio “ 7.5 logig Umax 


( 11 ) 


which gives 


1.8197 10 




23 


cm 


-31 




*4 


( 12 ) 


This formula is slightly different from the one proposed by Griem [37] (4 logj^g Z* instead of 4.5 logj^g 


in Eq. 0), who considered the electronic microfield instead of the ionic microfield for Eq. ( |TT| ). 

If the population Q of the Rydberg super-shell a is equal to zero, the opacity calculation involves only 
configurations whose contribution can be potentially detailed and therefore consists only of UTAs, SOSAs 
or DLAs. Indeed, a one-electron jump from an orbital (ni) gives rise to a transition array which can be 
detailed or not. In order to decide whether a line-by-line calculation is necessary and to determine the 
validity of statistical methods, it can be interesting to quantify the porosity (localized absence of lines) 
of transition arrays |38[ I39j . In SCO-RCG, the most limiting parameter is the maximum number of lines 
per transition array allowed by the user. There are also some computational constraints, mostly due to 
the sizes of arrays in the RCG routine: the number of lines must not exceed 800,000, the size of a J-block 
(number of levels having the same value of total atomic angular momentum J) of the Hamiltonian must 
not exceed 4,000 and we can not handle configurations with more than two electrons (or more than two 
vacancies) in any orbital with £> A (due to the tabulated coefficients of fractional parentage |4D1141j ~l. 
Optionally, to speed up calculations, it is also possible to resort to a criterion in order to check the 
coalescence of the transition array. For instance, we may define the parameter 



Wph 


(13) 


where Wph is the largest physical broadening (natural, Doppler, Stark, collisions, autoionization, etc.). In 
SGO-RCG, all the lines in a given transition array have the same physical broadenings and the average 
distance (in energy) d between two lines in the transition arrays is evaluated as 


d = 


V^UTA 

z 


(14) 


UUTA being the UTA variance [22j and L the number of lines, which can be evaluated exactly by efficient 
recursive techniques or estimated by an approximate formula |42] . If ^ is smaller than a particular value 
e with e << 1, the contribution of the transition array to opacity is calculated statistically within the 
UTA or SOSA formalisms. 

In one option of the code, if Q ^ Ij ^my transition starting from this super-configuration is treated 
in the STA formalism. However, such statistical modeling of Rydberg spectators may be irrelevant for 
some well-resolved transition arrays, as we will see in the next section. 


4 Test case: F-like 2p—> 4d transitions of iron. Is the statistical 
modeling of Rydberg spectators relevant? 

In the case of an iron plasma at T=182 eV and Nf.=3.1 10^^ cm“^ (conditions of a recent experiment 
performed by J. Bailey |5])|^ the statistical modeling of satellites [Be] 2p®(5s • • • 8d)^ —>• 2p^4d^(5s • • • 

^In the present paper, we will not show any experimental spectrum nor discuss any detailed aspects of the experiment 
(for further explanations see Ref. [21). 
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Figure 7: (Color online) DLA computation of the opacity of iron at r=182 eV and iVe=3.1 10^^ cm ^ 
without spectators and with spectators treated in a statistical way as a Rydberg super-shell (5s • • • 8d)^. 


Spectator n£ 

Number of lines 

P 

860 

d 

1506 

f 

1707 

g 

1724 


Table 1: Number of lines of transition array 2p®n£ —)■ 2p^ 4d^n£^, i varying from 1 to 4. 


8d)^ fill significantly the gaps around 11.4 and 11.5 A (see Figs. [^and|^ . To check whether it is relevant 
or not, we have studied the effect of each spectator electron separately (from 5s to 8d) on the spectrum. 


4.1 Study of the impact of spectators individually 

When the spectator electron moves from an s to an f orbital, the number of lines increases from 860 to 
1724 (see table f^. The histogram of the amplitudes (square root of the line strength) of the lines for 
2p®nd^ —>■ 2p^4d^nd^ with n=5 to 8 (see Fig. shows that the number of weak lines increases when 
the spectator moves far away from the nucleus (see table . The detailed calculation of satellites differs 
then from resonance lines only by a small shift and broadening (see Figs. 10 and 11). The shift of the 
center of gravity of transition array 2p^ —2p‘^4d^ due to the adding of a spectator electron is shown in 
Fig. It decreases globally as a function of n with a non-monotonic behaviour inside a given shell: a 


slight increase for s—^-p—:>d and then a decrease for I >3. We can see (solid line of Fig. 12) that the main 
contribution to the shift stems from the direct Slater integrals, for which the respective variations 
are displayed in Fig. 13 It is worth noting that the quantity {ni,2p) increases for .^ = 0,1,2 and 


decreases for i >3 inside a particular shell, while F^^\ni,4,d) is an increasing function of i inside a shell. 
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Figure 8: (Color online) Detail of Fig. [^in the range of transitions F-like [Be] 2p® —2p‘^ 4d^ 



Figure 9: (Color online) Histogram of the amplitude (square root of the line strength) for 2p®nd^ —> 
2p"*^4d^nd^ with n=5 to 8. The calculation was performed with Cowan’s code (RCN/RCN2/RCG). 
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Figure 10: (Color online) Oscillator strength of 2p® —>■ 2p‘* 4d^ with a spectator electron nd^, n=5 to 8. 
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Figure 11: (Color online) Opacity of transition array 2p® —>■ 2p'^ 4d^ without any spectator and with a 
spectator electron (from 5s to 8d). An instrumental resolution of 0.005 A was chosen for clarity 
reasons. 
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Figure 12: (Color online) Shift of the center of gravity (first-order strength-weighted moment) of tran¬ 
sition array 2p® —>■ 2p‘* 4d^ due to a spectator electron n(.^ (from 5s to 8d) compared to the difference 
(ni,4:d) — (n£,2p). The calculations were performed with Cowan’s code (RCN/RCN2/RCG) 

and checked with formula (21). 
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Figure 13: (Color online) Comparison of Slater integrals 2p) and 4d) for n£ varying from 

5s to 8d. The calculation was performed with Cowan’s code (RCN/RCN2/RCG). 
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Spectator i 

Number of lines with / < /max/1000 

5p 

201 

6p 

254 

7p 

286 

8p 

310 

5d 

342 

6d 

468 

7d 

569 

8d 

652 

5f 

530 

6f 

641 

7f 

745 

8f 

820 


Table 2: Number of lines of transition array 2-p^ni —>■ 2p"*^ 4d^n£^ for which the oscillator strength / is 
smaller than its maximum value divided by 1000. 


4.2 Failure of the statistical modeling of Rydberg spectators 

We have first calculated the contribution to the total opacity of F-like lines of ls^2s^2p® —)■ ls^2s^2p^4d^ 
(dashed line of Fig. 14). We have then calculated the contribution arising from Ne-like lines of the type 
ls^2s^2p®n£ —)■ ls^2?2p^4d^n.^ with ni ranging from 5s to the last bound orbital 8d. 

In order to investigate the contribution of satellites, we proceed as follows. Considering Ci = Co (niC)^ 
with C'o=ls^2s^2p® and riiii ranging from 5s to 8d, we compare, using Cowan’s code (RCN/RCN2/RCG), 
the pure DLA contribution to opacity {i.e. neglecting the contribution of the Rydberg electrons) of 
Co —t CQ=ls^2s^2p^4d^, evaluated as 


t{hh') = 


1 JV Tre^h 1 
47ren A me Z 


X! 9aje faJ^a'J'lpaJ^a'J'ihu) 

OcJ GCq 


with Qo=Q, to the “DLA-averaged” opacity of Ci —>■ (7'=ls^2s^2p'^4d(ni£i)^) calculated as 


i{hiy) = 


1 Af Tve^h 1 
47reo A me Z 
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9aje 


0(E^j 


(15) 


(16) 


Cg{Ci,-" ,Cjv} aJeC 


with Q = 10 and UiC ranging from 5s to 8d. Z is the partition function. As can be seen in Fig. 
the total contribution of satellites (solid line), calculated by Eq. (15), is very similar to the resonance 


transition (dashed curve) calculated by Eq. (16), except that it is slightly shifted and broadened. These 


spectator electrons are very weakly bound to the ion, so that they do not interact much with the core 
electrons and perturb the transition very weakly. Therefore, the “DLA-averaged” spectrum is very close 
to the resonance spectrum, i.e., the pure DLA calculation without any spectators. 

We then repeated the calculation by using the statistical approach to model the effect of the spec¬ 
tator. All orbitals from 5s to 8d are now gathered in a single super-shell and the resulting transition 
array ls^2s^2p® (5s • • • 8d)^ —?> ls^2s^2p^4d^ (5s • • • 8d)^ (see Fig. 15) is now represented by only two 
unresolved relativistic substructures 2 pi /2 4d3/2, 2 p 3/2 —t 4d5/2, while 2 p 3/2 —>■ 4d3/2 is not visible. 

Therefore, it can not reproduce the porosity of the transition arrays. Moreover and unfortunately, in that 
case, the maxima of the statistical spectral features correspond to minima of the detailed spectrum. In 
other words, the statistical calculation fills the gap between the lines more efficiently than the detailed 
one. 
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Figure 14: (Color online) Opacity of transition array 2p^ —>■ 2p‘* 4d^ without any spectator (DLA) 
and with a spectator electron (from 5s to 8d). “DLA-averaged” means that the spectrum was 


averaged over the different spectra corresponding to a different spectator electron, as described in Sec. 4.2 


The calculations were carried out with Cowan’s suite of codes (RCN/RCN2/RCG) and an instrumental 
resolution of 0.01 A was chosen for clarity reasons. 
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5 New approach 

Following the PRTA approach 1131231111 mentioned in the introduction, we suggest the replacement of 
the super-transition array S ^ S' 

f S = K2LS(3s)2(3p)6(3d)i0(4s)2(4p)3(4d)4(4f)6(5s5p • • • 8s8p8d)Q 
{ S' = K2L8(3s)2(3p)®(3d)9(4s)2(4p)3(4d)4(4f)7(5s5p • • • 8s8p8d)Q 

by the “reduced” transition array C ^ C' 

f C = K2L8(3s)2(3p)®(3d)i0(4s)2(4p)3(4d)4(4f)6 
{ C' = K2L8(3s)2(3p)6(3d)®(4s)2(4p)3(4d)‘‘(4f)7 

and under the assumption that the Slater integrals are the same for the initial and final (super-)configurations 
S and S' , the variance is additive with respect to the separation between active and passive subshells: 

v's^S' = [(5s5p5d5f5g • • • 8s8p8d)^] . (17) 

Therefore, the effect of spectators from 5s to 8d is included as an additional Gaussian broadening of 
variance Aspect [(5s5p5d5f5g • • • 8s8p8d)^] of the detailed lines of C —>■ C". 


5.1 Contribution of a spectator electron 


ance of transition array 


<>N' 

/3 


-)■ 




to the statistical shift and vari- 

7/N' + 1 

■/3 


The energy of a configuration C reads 


= + (iVy - (18) 

s s,s' 

where Ns is the electron population of subshell s and Ss,s' the Kronecker symbol. The one-body integrals 
are defined as 


Is — is + 


eV{r) — 


Ze 


2 1 


y1{r)dr 


47reor 

where j/s, radial part of the wavefunction multiplied by r, obeys Schrodinger equation 


hHji + 

2m dr'^ 2mr‘^ 


and the two-body interactions are given by 


1 ) 


— eV{r) 


ys{r) = isysir), 


(19) 


( 20 ) 


^'>'3 — s- 

9i - di- 


*J) - ^ XI ( 0 0 o) 


( 21 ) 


being the configuration-interaction radial integral of order k represents the 

direct Slater integral and R^^^ij,ji) = the exchange Slater integral). In the general case, the 

statistical shift due to Np spectators in subshell tp reads 




/C' - /C -p 

p p ' pp pen. 




-c 


pk . 


k^p 


+ -NpiNp-l) 


yC" _ yC 

* pp pp 


When Slater integrals are assumed to be the same for C and C, the shift becomes 


( 22 ) 
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where 


Np (Vpp — Vpa) , 


( 23 ) 


^/3 ^pa. 


F^°'> {nplp^nplp) - {riaia^npip) 

0 y G^^Hn^e^,npip) 


C £a k 

0 0 


<.p 

0 


(n„4 


I ^p-^p) 


(24) 


In practice, in SCO-RCG, this shift is not required, since each super-configuration has its own self- 
consistent potential and its own set of wavefunctions. 

For the contribution of spectator subshell £p^ to the variance, the general formula is more complicated, 
and can be found in Ref. [43]. When Slater integrals are assumed to be the same for C and C", the 
formula reduces to a more compact form and reads 


Ffp {iip + 2 Np) . (gg) 

{Up + 1) p ’ 


(25) 


where the two-electron variance = A {£p£a £p£p) is given in table II of Ref. |22| (formulas El to 

E8 and E7, p. 2430). 


5.2 Contribution of a spectator super-shell to the statistical shift and 
variance 

Let us consider a super-shell (for instance, in the preceding case, a = (5s 5p 5d 5f 5g • • • 8s 8p 8d). If 
neither a nor /3 do belong to cr, then the contribution of the “spectator” super-shell a to the center of 
gravity reads 


(26) 

n—1 ^ 

with 

cl,„,{D) = -J2{-Xsr9sDf, (27) 

s^cr 

where and = VsjS — Vsa- The usual STA partition functions Uq 

can be calculated by recursion relations (Hj . The variance to be included in the convolution of the DLA 
spectrum reads 


with 


»;r[ 


=i: 


Gq-n 

Uq 


Vn - 




n—1 

?7n = E 4>m{D)(l)n-m{D) + 71 {(pn {D"^) + On) , 
m—1 


where 4>m{D) is dehned by Eq. (27), 


(28) 


( 29 ) 
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Figure 15: (Color online) Opacity of transition array 2p^ —2p'^ 4d^ calculated by SCO-RCG without 
any spectator (DLA), with a spectator electron (from 5s to 8d) following the procedure described in 


Sec. 4.2 and using a statistical modeling. 




af3\ 


(30) 


SG<T 


and 


On = -y^{-Xsrgs 


sGcr 


(31) 


If a € cr or /3 e cr, then a statistical calculation is performed. We can see that the new modeling fills 


the gap much less than the previous approach (see Figs. 15 and 16). As expected, the structure is now 
close to the resonance one and slightly shifted. We can also observe that the opacity levels of the three 
substructures are different than the ones of the resonance. 

The algorithm of the SCO-RCG code including the modeling of Rydberg super-shell is described 
in Fig. 


5.3 SWAP approximation 

In the SWAP (Statistical Weight Approximation) limit, i.e. when temperature tends to infinity, one has 
Xg = 1 and ~ ^ Q ^ where G = 9s is the total degeneracy of super-shell a. After simple 

algebraic manipulation, one finds 




ia/3 


G 


for the shift of the detailed transition array a —)■ /? and 


(32) 
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Figure 16: (Color online) Opacity of transition array 2p^ —2p'^ 4d^ calculated by SCO-RCG without 
any spectator and within the previous and new modelings of Rydberg spectators. 


spect 


[Q] = 


Q{G-Q) 




q:/3 




(33) 


a/3 1) I J ' G{G-1) 

for the contribution to the variance. It can be checked in table that even at T=182 eV, the SWAP 
approximation for the correction to the shift and broadening due to Rydberg spectators is rather close 
to their exact values. 


5.4 Reducing the statistical part 

As explained in section in SCO-RCG, the total photo-excitation opacity consists in two parts: a 
detailed part, and a statistical part (which is expected to be much smaller than the detailed one). We 
can see for the case of Da Silva experiment [IS] at T=25 eV and p=0.008 g.cm“^ that the statistical 
part of the calculation has been significantly reduced within the new modeling of Rydberg spectators (see 
Figs. 18 and 19). The main An=0 (3-3) structure around hv=7{) eV is now closer to the experimental 
spectrum, although there are still some discrepancies (see Fig. 20). 
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Figure 17: (Color online) Algorithm of the SCO-RCG code. 
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Figure 18: (Color online) Detailed part, statistical part and total photo-excitation opacity calculated by 
SCO-RCG within the previous modeling of Rydberg spectators in the conditions of Da Silva’s experiment 
(r=25 eV and p=0.008 g/cm^) [45]. 
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Figure 19: (Color online) Detailed part, statistical part and total photo-excitation opacity calculated by 
SCO-RCG within the new modeling of Rydberg spectators in the conditions of Da Silva’s experiment 
(r=25 eV and p=0.008 g/cm^). 
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Transition 

Total variance (eV^) 

Contribution of (eV^) 

in SWAP (eV2) 

2s —> 3p 

0.3132 

6.42 10-4 

5.28 10-4 

2s -> 4p 

0.3518 

2.64 10-3 

2.29 10-3 

2p —>• 3s 

0.1312 

6.25 10-4 

5.15 10-4 

2p — y 3d- 

0.1129 

6.22 10-^ 

5.05 10-^ 

2p —)• 4s 

0.1530 

2.34 10-3 

2.04 10-3 

2p ^ 4d 

0.1419 

2.61 10-3 

2.24 10-3 


Table 3: Super-configuration (ls)^(2s)^(2p)^(j^ of iron at T =182 eV and 7Ve=3.1 10^^ cm 
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Figure 20: (Color online) Comparison between the experimental spectrum measured by Da Silva et al. 
and transmissions calculated by SCO-RCG within the previous and new modelings of Rydberg spectators 
in the conditions of Da Silva’s experiment (r=25 eV and yO=0.008 g/cm^) [45] , 


5.5 Impact on the Rosseland mean 

The new modeling may have a significant impact on the Rosseland mean, which is very sensitive to 
hollows and windows in the spectrum, and therefore to the line wings. Tableindicates the contribution 
of spectral band [1.5 - 2 keV] to Rosseland opacity, i.e. 


Kb. [hvi - = 


chu^ dB^(T) 
Ihvi dT 


dv 


rhu2 dB„{T) 1 ■ 

J hi^i dT K{hv) 


(34) 


with hvi=l.b keV and hi'2=‘2 keV, for a bromine plasma at T=200 eV and p=0.01 g.cm“^ calculated 
by SCO-RCG within three options: pure DLA (no Rydberg spectators), previous and new modeling of 
Rydberg spectators (Figs, [^and 22). It is clear that omitting the Rydberg electrons is not acceptable 
and that the differences between the values corresponding to the previous and new approaches are mostly 
due do the fact that the new model, which can be viewed as a perturbation to the pure DLA computation, 
fills the gaps between lines in a much less pronounced way than the previous statistical approach. 


5.6 Remark concerning autoionizing states 

Some autoionizing states [Ml HZ] may lead to an additional broadening of the absorption structures (not 
taken into account in SCO-RCG). For instance, the super-configuration of ion Fe^®’*' lies above 

the fundamental super-configuration of ion Fe^^+ (see Fig. 23). Therefore, considering only one 

electron in the Rydberg super-shell (which, in the present case, starts at orbital 5s), this effect is not 
important. However, for two or more electrons in the Rydberg super-shell, this broadening mechanism 
must be taken into account. Work is in progress to include it in SCO-RCG. 
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Energy (eV) 


Figure 21: (Color online) Bromine opacity calculated by SCO-RCG: pure DLA (no Rydberg spectators), 
previous and new modeling of Rydberg spectators. 



Figure 22: (Color online) Detail of Fig. 21 


28 



















































































KiLSQi 



K^L'^O^ 






Ground levels of Fe^^"^ 


K^LS 

Ground levels of Fe^^"^ 


Figure 23: (Color online) Autoionizing complexes in the case of iron. 
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Transition 

DLA 

(no Ryd. spec.) 

DLA 

+ stat. modeling 
of Ryd. spec. 

DLA 

+ new modeling 
of Ryd. spec. 

kr[1.5 — 2 keV] (cm^/g) 

901.27 

1227.52 

1129.42 


Table 4: Contribution of the frequency band [1.5 - 2 keV] to the Rosseland mean with the previous 
statistical model and the new approach compared to the calculation without DLA spectators. The 
Rydberg super-shell contains subshells 5s to 15t. 


5.7 Prescription for the inclusion of Rydberg spectators in opacity codes 

Let us consider the super-conhguration 

f SQ = CaQ 

\ C = ^ ’ 

a representing the Rydberg super-shell. One could consider convolving the detailed transition array 
C —>■ C" by a linear combination of profiles 

^ (hv - SEQ,a^^^' -bug) (36) 

where the shift reads 


the correction to the variance 


Q 

SEq = J2 


n—1 


Ug-nig) 

Uq{9) 


K{D), 


and 


is the Voigt function. 


Q 

H 

n—1 


UQ-n{g) 

Uqig) 


{SEgf 


H 


x.a.u) = 


rV2 


TTV J-c 


+ {x — J/)^ 


dy 


(37) 


(38) 


(39) 


6 Conclusion 

Spectator electrons in high-lying states may be important in X-ray absorption or emission spectra in LTE 
or NLTE hot plasmas. We proposed a new model, based on the additivity of the variances from active and 
spectator electrons in the STA theory. This approach, inspired from the Partially-Resolved-Transition- 
Array model, consists in a reduced detailed-line-accounting calculation, omitting these Rydberg electrons, 
the effect of the latter being included through an additional shift and broadening of the lines, expressed 
in terms of canonical partition functions. The resulting method can be used in any DLA opacity code. 


7 Appendix: The formula of Goldberg and Rozsnyai 

In their work Goldberg and Rozsnyai neglected the effect of terms as well as explicit electron-electron 
interactions in the Boltzmann factors. Their expression was obtained from the formula 
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^dielectronic — 


(40) 


E [(K) - (^p)'] > 

p Rydberg 
a,P 


where the average value of a function {A with N = (iVi, iV 2 , N^, - ■ ■ , N^^), is obtained as 

91 92 Qu, ( u 


(A (n)) 


E E- E n 

Ari=0W2=0 Af„=0 U=1 


5i 

iV,; 


X 


Ni 


N^{gp-Np)A(N) 


E E •E n 

Ni=0N2=Q N^=0 U=1 


9i 


(41) 


X 


Ni 


Na {gp - Np) 


with Xi = e This is a grand-canonical evaluation of the averages, i.e. without any assumption 


concerning the number of electrons. Applying formula (41) to A and A = X^ and using 

elementary binomial relations, one finds 


^dielectronic — 


E {Vpp-Vp^fN^(^ 

p Rydberg 

P ^ 06, ^ 


where Np = QpXpj (1 -|- Xp) is the average population of orbital p. 



(42) 
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